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ABSTRACT 


The  geometry  of  marine  propellers  is  described  in  this  report.  The  types 
of  data  used  for  blade  section  description  are  covered  and  methods  for  generating 
blade  section  leading  and  trailing  edge  forms  are  derived.  A  method  for  generating 
the  blade  tip  from  minimal  data  is  given.  The  computer  program  BLADE  developed 
for  generating  3  dimensional  propeller  blade  coordinates  is  described  together  with 
examples  of  its  graphic  output. 


SoTOTiaire 


Ce  rapport  decrit  la  geometrie  des  helices  marines.  Les  auteurs 
traitent  des  types  de  donnees  utilisees  pour  decrire  le  prof 11  de  la  lame  et 
le  developpent  des  methodes  de  production  des  formes  des  bords  d'attaque  et  de 
fuite  des  lames.  De  plus,  une  methods  de  production  de  1* incidence  de  la  lame 
a  partir  d'une  tres  faible  quantite  de  donnees  est  presentee.  Les  auteurs 
decrivent  aussi  le  logiciel  BLADE,  developpe  pour  le  calcul  des  coordonnecs 
tridimensionnelles  de  lame  d'helice  et  donnent  des  exemples  de  sorties 
graphiques . 
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NOTATION 


coefficient  of  cubic  spline,  basic  form  offset  at  point  i 
coefficient  of  cubic  spline,  slope  at  end  point  i 
coefficient  of  cubic  spline 

c  blade  section  chord  station  as  fraction  of  chord 
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C.  tabulated  section  design  lift  coefficient  (NACA  section  data)  per  unit 
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lift  coefficient  ratio 
r 

FR  fillet  radius 

LER  Leading  edge  radius  (NACA  section  data) 

th 

pitch  of  blade  at  the  i  section 

PMT  point  of  maximum  section  thickness  as  fraction  of  chord 
PT  Points  of  tangency  at  leading  edge 

RF  fraction  of  full  blade  radius  for  a  specified  section 

Rj  blade  section  radius  for  l^^  section 

RR  full  blade  radius  from  center  of  rotation  to  blade  tip 

RK  blade  rake 
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r^  radius  of  blade  tip 
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SK 
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max 


maximum  section  thickness 

thickness  of  the  blade  at  the  intersection  with  hub 

offset  distance  from  section  chord  to  back  surface 
offset  distance  from  section  chord  to  face  surface 
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trailing  edge  radius 
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chordline  distance  to  where  a  normal  to  the  meanline  meets  the  face 
contoxir  of  a  section 
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contour  of  a  section 
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X  coordinate  for  the  center  of  the  leading  edge  radius 

global  X  coordinate  for  blade  back  surface 

global  cartesian  coordinates  for  measurement  of  blade  back  surface 
global  cartesian  coordinates  for  measurement  of  blade  front  surface 
basic  form  offset  as  percent  of  chord 
mean  line  offset  at  x  modified  by  lift  coefficient 

mean  line  offset  at  x  as  percent  of  chord 
mean  line  offset  with  lift  coefficient  of  1 
maximum  value  of  y 
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mean  line  offset  as  fraction  of  chord 
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1.  INTRODUCTION 


The  Defence  Research  Establishment  Atlantic  (DREA)  is  conducting 
research  on  the  hydrodynamics  and  strength  of  marine  propellers.  Because  a  marine 
propeller  is  a  complex  three-dimensional  form,  the  shape  must  be  described 
adequately  for  strength  determination  or  manufacture.  The  approach  generally 
taken  is  to  define  the  geometry  of  a  single  blade  in  detail  including  the  intersection 
with  the  hub. 


Many  different  blade  shapes  have  been  developed  (Figure  1)  for  reasons 
such  as  the  reduction  of  unsteady  hydrodynamic  loading  on  the  blades  when  the 
propeller  is  operating  in  a  highly  non-uniform  wake  behind  the  ship's  hull.  Thus,  the 
method  chosen  for  geometry  definition  must  be  able  to  take  into  account  a  wide 
variety  of  blade  forms. 

Propeller  blade  geometry  is  most  often  supplied  as  two-dimensional  data 
in  the  form  of  blade  sections  from  which  three-dimensional  spatial  coordinates  must 
be  generated.  Specification  of  the  blade  shape  can  be  complicated  by  the  fact  that 
no  one  universal  format  and  axis  system  has  been  established  by  blade  designers  for 
presenting  the  geometry  data.  Therefore,  the  method  developed  to  generate  the 
spatial  coordinates  required  for  structural  analysis,  manufacture  and  measurement 
must  be  easily  adapted  to  the  various  forms  in  which  the  data  are  presented. 

This  report  describes  in  detail  the  geometry  of  propellers  and  an  exact 
method  for  converting  the  two-dimensional  data  into  three-dimensional 
coordinates.  The  method  can  handle  blade  geometry  in  the  form  of  NACA  section 
thickness  and  camberline  data  or  in  the  form  of  individual  section  coordinates.  In 
regions  such  as  the  leading  edge  and  the  blade  tip,  where  coordinate  information  is 
frequently  insufficient,  the  procedure  has  been  further  developed  to  fit  acceptable 
faired  curves. 

To  reduce  the  effort  required  to  produce  propeller  geometry  data,  a 
computer  program  has  been  developed  based  on  the  method  of  geometry  definition 
described  in  this  report.  It  can  be  \ised  to  generate  three-dimensional  blade 
coordinates  and  to  display  graphically  the  geometry  data  and  the  propeller  blade 
shape  according  to  naval  architectural  convention. 


2.  PROPELLER  GEOMETRY 
2.1  Blade  Description 

Marine  propellers  are  best  described  by  defining  their  parts.  The  four- 
bladed  variable  pitch  propeller  shown  in  Figure  2  illustrates  the  ma)or  propeller 
parts.  Variable  pitch  propellers  can  rotate  each  blade  around  its  spindle  axis; 
otherwise,  they  are  similar  in  geometry  to  fixed  pitch  propellers  except  at  the  blade 
hub  intersection.  This  intersection  portion  cf  a  variable  pitch  propeller  is  referred 
to  as  the  blade  palm,  which  is  a  circular  base  centered  on  the  spindle  axis.  The  palm 
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frequently,  though  not  always,  forms  a  portion  of  the  spherical  shaped  hub.  A 
sphere  is  used  to  ensme  a  good  hydrodynamic  fairing  at  the  blade  hub  intersection 
regardless  of  the  angle  at  which  the  blade  is  set. 

Fixed  pitched  propeller  blades,  on  the  other  hand,  (Figure  3)  intersect 
with  a  hub  that  can  vary  in  shape  from  spherical  to  cylindrical  and  from  barrel  to 
conical.  In  this  case,  the  blades  are  cast  integral  with  the  hub  or  they  are  bolted  to 
it. 


A  propeller  blade  has  two  main  hydrodynamic  surfaces.  The  surface  of 
the  blade  which  faces  aft  is  referred  to  as  the  face  of  the  blade.  The  surface  which 
faces  forward  is  called  the  back  of  the  blade.  A  propeller  which  normally  rotates 
clockwise  when  viewed  from  astern  is  called  right  handed  and  the  leading  edge  is 
that  edge  to  the  right  formed  by  the  intersection  of  the  face  and  back  surfaces  of  an 
upward  oriented  blade.  The  trailing  edge  Is  the  edge  to  the  left.  The  reverse  is  true 
for  counterclockwise  rotating  propellers  generally  referred  to  as  left  handed 
propellers.  The  convention  used  in  this  report  assumes  a  right  handed  propeller; 
therefore,  the  illustrations  of  the  geometry  are  shown  with  the  leading  edge  to  the 
right.  The  exceptions  are  in  the  derivation  of  section  coordinates  where,  for  more 
mathematical  convenience,  the  sections  are  shown  with  the  leading  edge  to  the  left 
in  accordance  with  airfoil  convention. 

The  tip  of  the  blade  joins  the  leading  edge  to  the  trailing  edge  at  the 
intersection  of  the  face  and  back  surfaces  at  the  maximum  radius  from  the  center 
of  rotation  from  which  the  blade  extends. 

To  make  a  smooth  transition  from  the  blade  to  the  hub  or  palm,  a  fillet  is 
used  (Figure  4).  The  fillet  is  most  easily  formed  by  a  series  of  single  radius  curves, 
the  radii  varying  between  a  maximum  near  mid  chord  and  a  minimum  at  the  leading 
and  trailing  edges.  Compound  curve"  are  often  xxsed  in  attempts  to  produce  constant 
stress  fillets.  These  are  described  in  more  detail  later  in  the  report. 

The  assembly  of  blades  forming  a  propeller  has  often  been  referred  to  as  a 
screw,  as  in  the  case  of  machine  screw,  because  they  are  similar  in  concept.  A 
propeller  blade  is  generated  from  a  series  of  fully  developed  sections  such  as  the  one 
shown  in  Figure  S.  Each  of  the  sections  forming  the  blade  is  located  at  a  specified 
radius  from  the  center  of  rotation.  Each  section  is  then  rotated  about  its 

vertical  or  spindle  axis  to  a  prescribed  angle  and  then  wrapped  or  curved  so  as  to 

fit  on  a  fictitious  cylinder  with  a  radius  equivalent  to  the  specified  radius  for  the 
section.  The  chord  line  of  the  section  thus  forms  part  of  a  helix  on  the  cylinder  in 
the  same  manner  as  a  screw  thread.  By  treating  each  of  the  sections  in  this  manner 
and  filling  and  blending  between  the  sections,  two  smooth  and  continuous  surfaces 
are  produced  forming  the  face  and  back  of  the  propeller  blade. 

2.2  Blade  Coordinate  Systems 

The  coordinate  system  for  a  propeller  blade  in  most  cases  is  specified  in 
the  form  of  right  handed  XYZ  cartesan  coordinates.  Unfortunately,  there  is  no 
uniform  standard  for  propeller  coordinate  systems,  as  shown  in  Figures  6a  to  6d. 

The  Defence  Research  Establishment  Atlantic  (DREA)  axis  system  (Figure  6a)  has 
the  Z  axis  positive  aft  coinciding  with  the  axis  of  propeller  shaft  rotation.  The  Y 
axis  is  the  vertical  reference  referred  to  as  the  pitch  or  spindle  axis,  about  which  a 
controllable  pitch  blade  is  rotated  to  change  pitch.  The  Y  axis  is  positive  up  and 
thus  the  X  axis  is  positive  to  starboard. 
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This  axis  system  was  chosen  because  it  maintains  the  XY  plane  as  the 
major  projection  plane  when  looking  forward  thus  placing  the  principal  transverse 
view  in  the  most  widely  used  XY  coordinate  plane.  Blade  cross  sections  are 
described  in  the  XZ  plane  which  is  close  to  the  convention  used  for  NACA  section 
data  where  X  is  the  chord  line  axis. 

2.3  Pitch 

As  described  previously,  a  propeller  blade  is  formed  from  fiiUy  developed 
sections  located  at  specified  radii  placed  so  that  their  chord  lines  are  at  angles 

to  the  horizontal  X  axis  in  the  XZ  plane  ^Fig^^^e  7).  This  angle  4>^  called  the  pitch 

angle,  follows  the  right  hand  convention  of  positive  angles  counterclockwise  about 
the  Y  or  pitch  axis.  A  view  of  the  section  on  the  XZ  plane  is  also  shown  in  Figure  8 
which  represents  a  surface  of  a  cylinder  of  radius  roiled  out  to  its  f\illy  developed 

length  2ttR^  .  When  this  surface  is  vinrapped  around  a  cylinder  of  radiiis  R^,  the  chord 

line  and  its  extension  form  a  true  helix. 

In  one  revolution  of  the  cylinder,  a  point  on  the  helix  travels  a  negative  Z 

distance  P.  =  2TrR.tan  0. 

1  1  1 

where  P.  =  pitch  of  section  i 

i 

R^  =  radius  of  the  cylinder  for  section  i 
0.  =  pitch  angle  of  section  i 

The  pitch  P^  is  the  same  axial  distance  a  screw,  of  the  same  cylindrical 
diameter  and  pitch  angle,  would  advance  in  one  revolution. 

The  propeller  designer  must  therefore  specify  the  pitch  Pj  or  the  pitch 
angle  0^  for  each  section  at  its  radius  R^.  He  must  also  specify  the  full  diameter  D 
of  the  propeller.  Frequently,  these  data  are  supplied  as  dimensionless  values  P^/D 
called  pitch  to  diameter  ratios  for  each  radius  R^.  The  pitch  angles,  if  not  given, 
can  be  extracted  from  this  information. 


=  tan  ^  — I -  (2.1) 

2TrRj/D 

A  pitch  diagram  (Fig\ire  9)  can  be  drawn  to  show  the  pitch  distribution 
with  respect  to  the  cylindrical  radii.  The  diagram  should  form  a  smooth  curve  when 
plotted  with  the  pitch  P^  as  the  abscissa  and  the  radii  as  the  ordinate,  thus  acting 

as  a  check  of  the  data  and  the  design. 
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Pitch  is  also  the  parameter  for  specifying  the  angle  to  which  a 
controllable  pitch  propeller  is  set.  Generally,  for  this  purpose,  the  pitch  angle  (t>.  of 

the  blade  section  at  the  radius  of  0.7  of  the  full  radius  is  used  as  a  reference 

angle,  and  the  setting  is  specified  as  a  P/D  ratio  where 


P  0.7  _  0.7Tr  tan  (J)  (2.2) 

^  - 


2.4  Skew 


Blades  of  simple  geometry  have  the  midchords  or,  in  some  cases,  points  of 
maximum  thickness  of  the  blade  sections  coincident  with  the  blade  pitch  or  spindle 
axis  (Figure  6).  Other  blades  of  more  sophisticated  geometry  have  their  sections 
displaced  along  their  chord  lines  from  the  spindle  axis.  This  displacement  is  called 
skew.  It  is  applied  to  reduce  the  unsteady  hydrodynamic  loading  on  the  blade  when 
the  propeller  is  operating  in  a  highly  non-uniform  wake  behind  the  ship's  hull.  Skew 
can  be  shown  on  the  fully  developed  unwrapped  view  of  the  blade  section  (Figure  10) 
as  a  linear  dimension  along  the  section  chord  line  from  the  intersection  with  the  Z 
axis.  The  skew  of  the  section  is  referenced  to  its  midchord  or  point  of  maximum 
thickness  (ie.  the  blade  section  reference  point).  Positive  skew  moves  the  section 
around  the  shaft  axis  in  the  flow  direction,  opposite  to  the  direction  of  shaft 
rotation. 


Skew  can  also  be  shown  as  a  true  linear  dimension  on  the  diagram  of  the 
expanded  view  of  the  blade  (Figure  11).  The  fully  developed  sections  are  laid  out 
flat  on  the  XY  plane  with  their  section  reference  points  displaced  in  the  direction  of 
the  negative  X  axis  a  distance  equal  to  the  skew  along  the  blade  reference  line. 
Positive  skew  is,  therefore,  measured  along  the  negative  X  axis  and  negative  skew 
along  the  positive  X  axis. 

Skew  is  sometimes  specified  in  angular  form  as  an  angle  a^,  called  the 
skew  angle,  projected  on  the  XY  plane,  through  which  the  wrapped  section  reference 
point  is  displaced  (Figure  12).  In  this  case,  the  linear  skew  dimension  SK^  for  the  i^^ 
section  can  be  obtained  by  multiplying  the  projected  skew  angle  (radians)  by  the 
radius  and  dividing  by  the  cosine  of  the  pitch  angle  ie. 


®i“i 


cos  (b 


I 


(2.3) 


2.5  Rake 

Rake  is  an  aft  or  forward  displacement  of  the  blade  section  reference 
point  from  the  XY  plane.  It  is  \xsed  to  increase  clearance  between  the  hull  and 
propeller.  As  shown  in  Figure  13,  a  blade  section  can  be  caused  to  have  rake  in 
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two  ways.  In  the  case  of  skewed  blades,  rake  is  induced  when  the  sections  are 
displaced  along  the  pitch  helix.  The  Z  distance  between  the  section  reference  point 
and  the  intersection  of  the  chord  and  the  YZ  plane  is  called  the  skew  induced  rake. 
The  conventional  method  to  introduce  applied  rake  (RK)  is  to  offset  each  blade 
section  an  additional  distance  in  the  Z  direction  from  the  X  axis.  The  rake  distance 
RK  is  from  the  X  axis  datum  (actually  the  XY  plane)  to  the  chord  line  YZ  plane 
intersection.  Offset  aft  produces  positive  rake  and  offset  forward  produces 
negative  rake. 

2.6  Warp 

When  a  blade  is  given  skew,  it  automatically  moves  the  blade  sections  aft 
along  the  helix  imposing  the  skew  induced  rake  on  the  blade  (Figure  13).  If  it  is 
desirable  to  maintain  the  blade  sections  in  the  plane  of  rotation,  then  the  blade 
sections  are  moved  back  by  imposing  rake  on  the  sections  equal  to  the  negative  of 
the  skew  induced  rake  (Figure  14a).  Warp  can  thus  be  defined  as  the  angular 
displacement  equal  to  the  projected  skew  angle  of  the  blade  sections  measured  from 
the  Y  axis  with  any  axial  displacement  such  as  skew  induced  rake  removed  (Figure 
14b).  The  linear  dimension  for  the  warp  of  section  i  (Figure  14a  and  14c)  on  the 
\anwrapped  view  of  the  blade  and  the  expanded  outline  is 

WK.  =  R,a.  (2.4) 

1  i  1 


where  =  warp  or  the  projected  skew  angle  of  the  i 


th 


section. 


2.7  Blade  Sections 


The  logical  process  for  determining  the  geometry  of  a  propeller  begins 
with  the  selection  of  the  blade  section.  The  designer  may  develop  his  ovni  or  choose 
from  a  variety  of  proven  sections  whose  hydrodynamic  performance  have  been 
characterized  and  recorded. 

A  blade  section  (Figixre  15)  is  initially  described  in  its  fully  developed  or 

expanded  form^.  In  this  form,  the  shape  is  most  often  described  by  offsets  TF  and 
TB  at  fractions  c  of  the  section  chord  line  length  S.  The  section  contour  formed  by 
the  upper  offsets  TB  measured  from  the  chord  line  or  datum  is  called  the  back  or 
suction  side  of  the  section.  The  section  contour  formed  by  the  lower  offsets  TF  is 
called  the  face  or  pressure  side.  The  chord  line  in  addition  to  forming  the  datum  for 
the  section  offsets  acts  as  the  line  by  which  the  section  is  angled  to  the  horizontal 
datum,  referred  to  as  the  pitch  angle.  Normally,  the  face  of  the  section  is  convex 
with  the  offsets  below  the  chord  line  which  are  considered  to  be  positive  values. 
When  the  offsets  are  all  on  or  above  the  chord  line  such  that  the  face  is  flat  or 
concave,  then  the  chordline  is  also  referred  to  as  the  face  pitch  line  (Figure  16). 
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2.7.1  Airfoil  Section 


If  airfoil  sections  are  chosen  for  marine  propeller  blades,  then  airfoil 
convention  is  xised  with  the  leading  edge  to  the  left,  and  the  offsets  are  given 
normal  to  a  camber  or  mean  line  rather  than  normal  to  the  chord  line  (Figtire  17). 
When  blade  section  data  are  presented  in  the  format  for  NACA  airfoil  sections,  the 
data  have  to  be  converted  to  offsets  from  the  chord  line.  The  basic  thickness  forms 
and  mean  line  data  for  a  series  of  sections  can  be  found  in  Reference  3. 


In  addition  to  the  basic  form  and  the  mean  line  data,  the  lift  coefficient 
and  the  leading  edge  and  trailing  edge  radii  are  also  supplied.  For  example,  a 
section  from  the  NACA  five  digit  series  may  be  specified  as  NACA  16-309  with  a 

mean  line  NACA  x  =  1.0  where  I  is  the  fraction  of  the  chord  uniformly  loaded  at  the 
ideal  angle  of  attack.  This  is  interpreted  as  a  16  series  basic  form  with  a  thickness 
to  chord  ratio  T/S  of  .09  or  nine  percent  as  Indicated  by  the  last  two  digits.  The 
first  digit  after  the  dash  is  the  lift  coefficient  Cj^  expressed  as  a  decimal.  The  lift 

coefficient  in  this  example  is  0.3.  Referring  to  Figure  17  and  page  382  of  Reference 
3,  the  mean  line  curve  is  plotted  from  the  data  table  or  from  an  equation  such  as  the 
following,  which  is  for  a  uniform  chordwise  pressure  distribution  at  zero  angle  of 
attack. 


It 

c  ytn  = 

-  ^ 

((1-  *m  )  In 

(l-*m)+*m  In  *m)  (2.5) 

dVc 

11 

4tt 

tin  (1 

S 

X, 

-  _L)  -  In 

S  S  S 

( ^)]  (2.6) 

dx. 

4^ 

S 

S 

where 


=  fixed  point  on  the  x  axis 


X  =  distance  along  chord  from  the  leading  edge 
m 


y  =  mean  line  ordinate  at  lift  coefficient  of  1.0 
'm 


x^  =  blade  section  chord  station  distance  as  percent  of  chord 


y^  s  mean  line  offset  ordinate  at  x  as  percent  of  chord 


y^  =  mean  line  ordinate  as  fraction  of  the  chord 


dy 

— s  slope  of  the  mean  line  at  x  _ 

dx  ^ 

c 

S  s  length  of  chord  line 


Cj  =  lift  coefficient 


A  mean  line  equation  which  may  be  more  s\iitable  for  marine  propellers  but  is  too 
complex  to  describe  here  can  be  found  on  page  74  of  Reference  3  as  equation  4.27. 
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The  basic  thickness  form  data,  page  309  Reference  3,  is  applied  to  the 

mean  line  ordinates  and  normal  to  the  mean  line  local  slope  dy  /dx  .  Other  design 

c  c 

lift  coefficients  can  be  accommodated  by  multiplying  the  tabulated  values  of  y  and 

c 

dy  /dx  by  the  ratio  of  the  desired  lift  coefficient  to  the  given  lift  coefficient, 
c  c 

The  normals  to  the  mean  line  are  then  transformed  to  offsets  from  the 
chord  line  for  the  upper  (back)  and  the  lower  (face)  surface  lines  of  the  section  by 
the  following  expressions.  The  expressions  are  based  on  a  local  axis  system  that 
agrees  with  airfoil  convention  (Figure  17).  The  expressions  for  TB  and  TF,  the  back 
and  face  offsets  are  introduced  at  this  point  as  the*e  symbols  were  used  previously 
for  non  airfoil  sections.  In  addition,  their  use  allows  the  algebraic  signs  for  the 
offsets  to  be  made  consistent  with  those  chosen  for  the  blade  axis  system  where 
offsets  above  the  chord  line  are  negative  and  below  the  chord  line  are  positive.  This 
approach  thus  allows  the  data  tables  from  Reference  3  to  be  used  directly  to  obtain 
TB  and  TF  values. 


© 

_  tan  ^  (Cj^  .  dy  /dx  ) 
r 

(2.7) 

r 

m 

^a 

(2.8) 

*m 

=  X  .s/100 
c 

^t 

=  y.S/100 

a  y^.  Cj  .S/100 

(2.9) 

=  X  -  y^  sin© 
m  t 

(2.10) 

=  y«  +  y  cos© 

'a  t 

(2.11) 

TB 

•-’'u 

(2.12) 

=  X  +  y^  sin© 
m  t 

(2.13) 

t- 

=  y^-y^  cos© 

(2.14) 

TF 

■ 

(2.15) 
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For  the  situation  where  the  data  tables  of  Reference  3  do  not  provide  the 
desu?d  thickness  to  chord  ratio  T  ,  the  basic  form  offsets  can  be  scaled  as  follows 


yt  = 


y 

^ax 


•S  T, 


(2.16) 


The  expressions  for  calculating  the  section  offsets  when  using  a  varying  Tj  ratio 
become 


11 

X 

T  «inA 

(2.17) 

m 

2y.  * 

^tmax 

Y  = 

V  4 

^t  T 

(2.18) 

U 

^  a 

\  = 

X  4 

^t  J 

(2.19) 

m 

2y*  ® 

tmax 

Y  = 

V  - 

y 

t  T 

(2.20) 

L 

2y»  * 

nmax 

2.7.2 

Interpolation  of  Section  Coordinates 

Interpolation  of  additional  section  coordinates  is  frequently  required.  One 
reason  is  that  the  upper  and  lower  offsets,  when  converted  to  the  chord  line  datum 
from  NACA  section  data,  do  not  always  originate  from  the  same  chord  locations  for 
each  section.  Though  in  this  form  section  coordinates  are  useful  for  graphically 
presenting  the  geometry  of  the  blade  in  spatial  X,  Y,  Z  cartesian  coordinates,  they 
are  not  suitable  for  stress  analysis,  particularly  when  using  the  fiiiite  element 
method.  It  is  desirable  therefore  to  generate  upper  and  lower  section  offsets  TB, 
and  TF  at  identical  chord  locations. 

A  computer  program  BLADGM  has  been  developed  for  generating  the 
desired  section  coordinates  and  is  described  in  Section  4  of  this  report.  The  program 
uses  cubic  spline  interpolation  which,  together  with  a  section  plotting  program,  is 
available  to  display  fully  developed  blade  sections  on  a  Tektronix  graphics  display 
terminal. 
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2.7.3  Leading  Edge 


The  leading  edge  radius  (LER)*  Figure  (18)i  is  generally  given  with  the 
basic  form  data.  The  radius  specified  in  the  basic  form  data  is  the  hypothenuse  of 
the  triangle  at  the  leading  edge  with  an  included  angle  6,  where  6  is  the  slope  of 
the  mean  line  at  an  X  value  of  O.S  percent  of  the  chords.  Thus,  the  coordinates  for 
the  center  of  the  leading  edge  radius  are 

X  =  LER  .  cos  [tan-1  (  .  *^^c  ))  (2.21) 

Yj  =  LER  .  sin  [tan‘^  (^4  .  ^  )]  (2.22) 

dXg  0.5% 

The  trailing  edge  radius  if  not  otherwise  specified  is  given  by  TER  =  .OIT 
where  T  is  the  maximum  thickness  of  the  basic  form. 

Often  insufficient  data  arc  supplied  to  accurately  describe  the  leading 
edge  region  of  a  propeller  blade  section.  There  are  a  number  of  methods  that  can  be 
used  for  generating  the  curve  through  the  leading  edge  from  the  data  supplied  for 
the  section.  One  such  method  depends  on  estimating  the  points  of  tangency  (PT)  of 
the  leading  edge  arc  (Figure  18)  and  the  curves  through  the  remaining  points  of  the 
back  and  face  contours.  This  method  is  awkward  to  program  and  depends  on  trial 
and  error  fitting  of  the  ctirves  using  a  spline  fit. 

2. 7. 3.1  Transformation  Method 

A  second  and  better  method  uses  the  same  basic  thickness  form  data  and 
the  leading  edge  radius  LER  to  fit  a  smooth  curve  through  the  data  points  and  the 
leading  edge.  A  cubic  spline  fit  was  chosen  for  this  method  as  it  is  compatible  with 
the  thickness  form  equations  given  on  page  117  of  Reference  3. 

The  spline  fitting  technique  is  based  on  passing  a  set  of  cubic  functions 

through  the  points  using  a  new  cubic  in  each  interval.  The  equation  for  the 

interval  between  points  x,  and  x,  ,  is* 

i  i  +  1 

y  =  A^  +  Bj  (x-Xj )  +  Cj  (x-Xj)*  •»  (x-x^)*  (2.23) 


The  cubic  must  fit  the  two  end  points  of  the  interval  x^  S  x  S  x^  ^  ^ 
(Figure  19a).  Therefore  at  x  »  x^ 

=  Aj  (2.23a) 

Simularly  at  x  =  x^^  ^ 

y,  ,  *  A,  4  B,  (x,  ,-x,)  4  C,  (x,  ,-x,)*  4  D,(x,  ,-x,)*  (2.23b) 

'Ul  i  i  Ul  i  i  i4l  i  i  i4l  i 
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An  additional  requirement  is  that  there  must  be  slope  compatibility  at  the 
end  points  of  adjoining  intervals.  For  the  interval  £  x  £  x^^^^ 


at  X  =  X.  =  B,  {2.23c) 

1  dx  i 

at  X  =  X,  ,  M  ^  B,  +  2  C,  (X,  ,  -  x.)  4  3D,(x,  x,)^  (2.23d) 

i4l  ^  i  1  Ul  i  i  i4l  i 

the  A  ,  B.,  C.,  D.  coefficients  can  be  obtained  from  the  data  supplied  provided  the 
1  1  1  i 

offsets  y.  and  the  slopes  can  be  found  at  the  end  points  of  the  intervals. 

Unfortunately,  at  the  leading  edge,  fFlgiire  19  (a))i  the  slope  is  infinite 

and  cannot  be  used  to  solve  for  the  coefficients.  This  problem  can  be  overcome  by  a 

2 

transformation  of  the  x  axis  in  the  interval  x  SxSx  of  x  =  C  •  The  derivation  of 

1  2 

the  transformation  is  given  in  Appendix  A. 

For  the  transformed  coordinate  system,  the  spline  equation  is 


y  =  A  4  B  (C~C  )  +  C  (C-C  4  D  (C-C  )*  (2.24) 

atC  =  C^  =  0  (2.24a) 

y  =  A  =  0  (2.24b) 

1  1 

B  (2.24c) 


dC 

where:  r  =  specified  leading  edge  radius  LER. 

atC  =  C 

2 

y  =B  C  +  C  C^4D  C’  =  A  (2.24d) 

^2  12  12  12  2 

^=B  4  2C  C  +3D  C’»2BC  (2.24e) 

1  12  12  22 
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As  can  be  seen  from  Figure  19(b),  the  transformed  curve  now  has  less  than 

infinite  slope  at  f  =  0  and  A  ,  B  ,  C  and  D  can  be  obtained  in  the  interval 

111  1 

Therefore,  since  C  =  vTc 

A  *  0  (2.25) 

1 

B  »  ■/Tr  (2.26) 


,  3A 

C  =  -L-  ( — 2-  -  2B  -  2B  ] 

^  Vx  VTT  2  2  1 

2  2 

1  2A 
D  =  J-  (  2B  VTC  -  — 2-  4  B  ] 

lx  3  2^  1 

2  2 

and  In  the  first  interval 

y  =  A4BV1c4Cx4D  (ymf 
11  11 

where  r  =  leading  edge  radius  LER 

X  =  end  point  of  interval  x  ^  x  S  x 
2  12 

A^  =  the  basic  form  offset  at  x^  generally  at  0.5  percent  of  chord 

B  =  the  slope  at  the  end  point  x 
2  2 

The  slope  is  obtained  from  a  cubic  spline  fit  of  the  offsets  y^  to  y  of 

the  basic  form  data  in  the  intervals  x  to  x  .  Additional  ordinate  values  are 

2  n 

thereby  obtained  in  this  first  Interval  from  equation  2.28a.  The  cubic  spline  fitting 

equation  defining  the  section  from  the  leading  edge  over  the  interval  x  to  x  has 

1  2 

been  programmed  and  incorporated  in  the  program  BLADGM. 

2. 7. 3. 2  Elliptic/Hyperbolic  Leading  Edge 

A  third  method  for  forming  the  leading  edge  is  to  use  an  ellipse  or 
hyperbola  (Figure  20).  In  this  method  the  leading  edge  is  formed  from  two 
elliptic/hyperbolic  portions  for  which  the  semi-axes  may  be  different;  one  for  the 
face  and  one  for  the  back.  The  ellipses  or  hyperbolae  (face  and  back)  must  have  the 
following  properties: 

(a)  Their  axes  must  be  aligned  with  the  chord  line, 

(b)  They  must  pass  through  the  leading  edge  point, 

(c)  They  must  pass  through  the  offset  y^  at  a  point  x^  nearest  to  the 
leading  edge,  and 


(2.27) 

(2.28) 

(2.28a) 


-11- 


(d)  The  slope  of  the  ellipse  or  hyperbola  at  the  offset  nearest  to  the 

leading  edge  is  the  same  as  the  slope  of  the  smooth  curve  at  that 

point. 

To  form  the  ellipse  or  hyperbola,  offsets  at  the  leading  edge,  which  are  a 
smooth  continuation  of  the  section  form  (Figure  20),  m\ist  be  given.  These  blade 
edge  offsets  are  used  together  with  the  other  available  offsets  to  obtain  the  slope  of 
the  section  form,  which  is  the  coefficient  of  the  cubic  spline  passed  through  the 

points  to  satisfy  (d).  If  the  blade  edge  offsets  are  not  available,  then  the  value  of 
leading  edge  radius  may  be  substituted. 

The  basic  equation  for  generating  the  elliptic/hyperbolic  leading  edge 
coordinates  is  developed  from  the  basic  equation  for  ellipses  or  hyperbolae, 


(a  -  ,  1  (2.29) 

2  .a 

a  b 

The  slope  at  the  offset  nearest 

bl 
2 

a 

B^  =  the  slope  of  the  spline  function  fitted  through  the  upper  or  lower 
blade  edge  point. 


the  leading  edge  is 


a-x 

( - i) 


(2.30) 


From  the  ordinate  and  slope  at  the  semi  axes  can  be  determined  as 


a  = 


X  (B  X  ♦ 
1X1 


2B  X  4  y 
11  1 


and 


(2.30a) 


b*=  y  (y  -  B  (a  -  X  )) 

^1  ■'i  1  1 


(2.30b) 


to  give  the  final  equation  for  generating  additional  coordinates  for  the  leading  edge 

in  the  interval  0  £  x,  5  x 
i  1 


^i 


2 

a 


(2.31) 
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2. 7.3.3  Leading  Edge  Extension  and  Final  Coordinates 

When  the  section  has  camber,  the  leading  edge  will  extend  beyond  the 

basic  form  chord  station  x  =  0  (Figures  18  and  21).  Fitting  the  leading  edge  to  an 

1 

already  cambered  section  will  introduce  more  complications  than  )ust  c\irve  fitting 
the  leading  edge  to  the  basic  form  before  camber  is  added.  If  not  already  given,  the 
offset  coordinates  have  to  be  shifted  temporarily  to  the  line  running  through  the 
extreme  leading  edge  point  where  there  is  a  leading  edge  radius  r.  The  chord  length 
must  be  increased  by  the  amount  it  extends  beyond  station  to  produce  the 

modified  chord  length  S  which  will  become  the  datum  for  all  the  chord  stations 

m 

and  section  offset.  Then  the  back  and  face  contours  are  curve  fitted  separately  for 
the  contour  of  the  nose.  The  interpolation  stations  for  identical  back  and  face 
positions  can  be  input  as  selected  values  or  generated  within  the  program  by  the 
following  equation  where  0  is  incremented  from  0®  to  180*  according  to  the  required 
number  of  stations. 


X  =  ^-OL  {  1  -  cos  0,  )  (2.32) 

i  2  i 

After  curve  fitting  the  nose,  the  vertical  coordinates  are  shifted  back  to  the  original 
chord  line  position. 


2.8  Expanded  View 

The  expanded  view  for  a  right  handed  propeller  (Figure  11)  has  been  used 
previously  to  illustrate  the  skew  of  the  blade  sections  of  a  propeller.  The  view  also 
shows  the  fully  developed  blade  sections  revolved  about  their  chord  lines  to  form 
projections  on  the  XY  plane.  The  sections  are  in  true  proportion  at  their  relative 
radial  positions  with  the  back  surfaces  of  the  sections  up  and  the  faces  down.  Their 
section  reference  points  are  aligned  with  the  Y  axis  or  a  specified  X  distance  from 
it  equal  to  the  skew,  if  it  is  present.  If  the  section  coordinates  have  been  specified 
from  the  leading  edge  as  fractions  of  the  chord  CFj,  then  they  are  transformed  to 

the  new  origin  as  follows  (Figure  22a).  The  chord  location  using  mid  chord  as 
datum  for  the  i^^  blade  section  is  given  by 

X^j  *  S^/2  -  (CFj  +  SKj)  (2.33) 

or  using  the  point  of  maximtun  thickness  as  datum  expressed  as  a  fraction  of  the 
chord  length 

Xjj  =  Sj.  PMT^  -  (CFj  Sj  ♦  SKj)  (2.33a) 
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where 


Xij  =  offset  station  for  section 

=  chord  length  for  i^^  section 

CF.  =  chord  fraction  for  offset 
) 


PMT^  =  point  of  maximum  thickness  as  new  datum  for  i  section  as  a 

fraction  of  the  chord 

SK.  =  skew  (see  Section  2.4) 

1 

A  further  complication  can  arise  if  the  chord  stations  CF^  are  specified  as 

fractions  of  the  dimensions  from  the  PMT  datum  to  the  leading  and  trailing  edges 
(Figure  22b).  For  this  case,  the  chord  stations  are  calciilated  in  two  steps.  From 
leading  edge  to  PMT 

X..  =  S,  .  PMT,  .  CF.  -  SK,  (2.33b) 

1)  i  i  )  i 

From  PMT  to  trailing  edge 

X,.  =  -(S,  (1-PMT,)  CF,  ♦  SK,)  (2.33c) 

i)  i  i  )  i 

where  "  fraction  of  chord  length  from  PMT  datum  to  leading  and  trailing 

edges 

The  sections  are  thus  arranged  in  the  proper  relationship  to  the  blade 
datum  which  is  the  Y  or  pitch  axis.  A  continuous  curved  line  drawn  through  the 
leading  and  trailing  edges  of  the  sections  forms  what  is  called  the  expanded  outline 
of  the  blade  projected  on  the  XY  plane  (Figure  11).  This  is  generally  chosen  as  the 
first  graphic  representation  of  the  propeller  blade  geometry  and  gives  the  chords, 
section  offsets,  section  radii  and  skew  in  true  proportion  in  accordance  with  the 
scale  of  the  drawing. 

2.9  Generating  the  Blade  Tip 

2.9.1  Tip  Outline 

The  geometric  details  of  the  propeller  blade  frequently  do  not  give 
sufficient  details  of  the  tip  outline  curve.  The  last  radius  and  section  specified  are 
generally  at  about  0.95  of  the  full  radius  of  the  propeller  (Figure  23).  It  is  necessary 
therefore  to  complete  the  blade  outline  at  the  tip.  This  is  not  a  difficult  job  for  the 
pattern  marker  or  the  draftsman  as  they  can  fit  a  curve  by  eye  using  drawing 
Instruments.  This  process  is  not  practical  for  a  computer  geometric  definition,  and 
another  method  must  be  used.  The  method  chosen  is  similar  in  most  respects  to  that 
used  for  ctirve  fitting  the  blade  section  leading  edge  as  described  in  2.7.3  of  this 
report.  In  this  case,  the  expanded  outline  is  used  except  that  skew  is  not  included. 
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The  outline  is  generated  from  the  end  points  of  the  chords.  The  chord  mid  points  or 
points  of  maximum  thickness  are  located  on  the  spindle  or  Y  axis.  Because  the  tip 
radius  is  not  specified  and  is  a  requirement  for  the  solution,  it  can  be  shown  (see 
Appendix  B)  that  based  on  a  parabolic  fit  and  by  using  the  leading  and  trailing  edge 
coordinates  of  the  last  given  section,  the  radius  of  curvature  at  the  tip  can  be  given 
by 


where 

x^  =  the  distance  from  the  leading  or  trailing  edge  to  the  Y  axis  at  the 
last  given  section  radius. 

y  =  distance  from  the  last  section  radius  R,  to  the  full  blade  tip  radius 
1  i 

RR. 

The  tip  outline  is  then  generated  using  the  cubic  spline  equation  described 
in  Section  2. 7. 3.1.  For  this  operation,  it  is  convenient  to  shift  the  origin  of  the  axis 
system  to  the  tip.  In  the  final  calculations,  the  tip  coordinates  are  transformed 
back  to  the  original  blade  axis  system  and  the  skew,  if  any,  is  added. 

2.9.2  Tip  Thickness 

Once  the  blade  outline  at  the  tip  has  been  formed,  then  the  blade 
thickness  from  the  last  specified  radius  to  the  tip  must  be  determined.  This  is 
accomplished  by  using  the  offsets  TB  and  TF  at  mid  chord  or  at  the  point  of 

maximum  thickness  for  each  of  the  specified  blade  sections.  The  offsets  when 
plotted  produce  a  section  which  is  a  good  spanwise  representation  of  the  blade 
section  shape  ]ust  as  the  horizontal  blade  sections  are  of  the  chordwise  shape 
(Figure  24).  From  these  data  the  section  at  the  blade  tip  is  generated  in  radial  or 
elliptic/hyperbolic  form  in  the  same  manner  as  the  leading  edge  of  the  chordwise 
sections  in  2. 7.3.1  and  2. 7. 3. 2.  If  the  radial  form  of  tip  profile  generation  is  chosen 
then  a  radius  for  the  tip  profile  must  be  specified.  If  an  elliptic/hyperbolic  profile  is 
chosen  then  a  knovni  or  estimated  thickness  for  the  semi-minor  axis  must  be  given 
for  generating  the  axes  and  thickness  offsets  using  the  equations  presented  in 
Section  2. 7.3.2.  In  the  case  where  the  center  of  the  blade  tip  (without  rake  and 
skew)  does  not  coincide  with  the  vertical  reference  axis  a  value  CT  for  the  offset 
may  be  given.  If  it  is  not  available  then  it  must  be  estimated.  (Figure  24). 

On  occasion,  the  computer  generated  drawing  of  the  blade  thickness 
diagram,  iwithout  rake,  may  show  the  blade  tip  profile  hooked  to  the  right  or  left. 
This  can  be  removed  by  experimenting  with  values  for  CT. 

2.9.3  Tip  Section  Data 

Blade  chordwise  section  data  for  the  additional  blade  radius  fractions  at 
the  tip  are  obtained  by  scaling  the  last  given  section  data  by  the  ratios  of  the  new 
and  given  chord  length  and  the  ratios  of  the  new  and  given  maximtun  section 
thicknesses.  Figure  25  shows  these  data  plotted  in  the  form  of  expanded  sections 
from  the  last  given  section  to  the  blade  tip.  Pitch,  skew,  and  rake  for  the  additional 
chord  lengths,  are  obtained  by  interpolating  the  original  data  provided  values  for 
these  parameters  are  given  for  the  tip  in  the  original  data. 
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2.10  Blade  Fillets 

The  transition  between  blade  and  hub  is  made  by  fitting  a  fillet  tangent  to 
the  blade  svirface  and  to  the  hub  or  palm.  A  variety  of  shapes  are  used  for  propeller 
hubs  and  palms.  For  variable  pitch  propellers,  the  palm  is  sometimes  of  spherical 
form  (Figure  3).  Using  the  same  spherical  shape  for  the  palm  guarantees  a  smooth 
uninterrupted  surface  at  the  blade  hub  regardless  of  blade  pitch  angle.  This  is  not 
aways  achievable  for  various  reasons,  in  which  case,  the  spherical  radius  for  a  palm 
can  be  larger  than  the  radius  of  the  hub. 

In  the  case  of  fixed  pitch  propellers,  the  forms  most  generally  used  for 
hubs  are  cylindrical  and  conical.  Distorted  versions  of  these  forms  are  called  bullet 
and  barrel  shaped  hubs  (Figxire  3). 

The  requirement  to  fit  a  fillet  for  smooth  transition  between  the  hub  and 
blade  for  such  a  variety  of  hub  forms  is  most  easily  accomplished  by  the  pattern 
maker  in  the  case  of  cast  propellers.  It  is  much  more  difficult  to  do  in  the  case  of 
blades  machined  on  an  N.C.  mill.  For  N.C.  milling,  all  surfaces  must  be 
mathematically  described  so  that  machine  tool  paths  can  be  generated. 

The  proportions  of  the  fillets  are  based  on  hydrodynamic  and  strength 
requirements,  the  idea  being  to  generate  a  fillet  of  such  proportions  as  to  give  a 
constant  stress  through  the  transition  while  meeting  hydrodynamic  and 
manufacturing  requirements. 

One  standard  for  a  constant  stress  fillet  proportions  in  common  use  is  a 
compound  radius  3Tj  and  T^  /3  where  T^  is  the  blade  thickness  at  its  Intersection 

with  the  hub  (Figure  26).  Unfortunately,  this  fillet  form  is  not  suitable  for  N.C. 
machining  because  it  is  discontinuous  in  the  second  derivative  at  the  tangency  points 
of  the  two  radii  and  the  tangency  points  with  the  blade  and  hub.  This  problem  has 
been  overcome  by  the  use  of  fillets  generated  by  a  continuous  function. 

2.11  Data  Verification 

Data  plots  of  the  blade  data  can  be  generated  by  the  computer  program 
BLADIN.  These  plots  allow  Inspection  of  the  supplied  tabular  data  in  graphic  form 
so  that  any  discontinuities  in  the  basic  data  can  be  detected.  Computer  generated 
plots  of  the  following  data  can  be  produced; 

(a)  Pitch  angle,  q>,  versus  blade  section  radius,  (Figure  27); 

(b)  Skew,  SK,  versxis  blade  section  radius,  (Figure  28); 

(c)  Rake,  RK,  versus  blade  section  radius,  R^  (Figure  29); 

(d)  Section  chord,  S,  versus  blade  section  radius,  R^  (Figure  30); 

(e)  Leading  edge  radius,  LER  versus  blade  section  radius,  R^  (Figure  31); 

(f)  Trailing  edge  radius,  TER,  versus  blade  section  radius,  R^  (Figure  32); 
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(g)  Maximum  camber,  y  ,  versus  blade  section  radius,  R,  (Figure  33); 

max  i 

(h)  Maximum  thickness,  T,  versus  blade  section  radius,  R^  (Figure  34); 

(i)  Individual  section  plot  (Figure  35). 

2.12  Generation  of  Spatial  Coordinates 

The  spatial  coordinates  for  the  blade  are  generated  by  transforming  the 
two  dimensional  coordinates  of  the  expanded  view  and  the  blade  sections  into  the 
three  dimensional  coordinates  of  the  blade  in  its  final  form  (Figure  36).  The  fully 
developed  sections  in  the  expanded  view  (Figure  11)  show  skew  as  an  offset  SK^  of 

th;  mid  chord  from  the  Y  axis.  As  stated  previously,  in  some  cases,  the  point  of 
maximum  thickness  PMT^  is  used  rather  than  the  mid  chord  point.  If  the  blade  has 

rake,  the  unwrapped  sections  on  the  XZ  plane  are  offset  a  distance  along  the  Z  axis, 
an  amount  equal  to  the  rake  RK^  (Figure  13).  Each  section  is  rotated  about  the 

point  of  intersection  of  the  chord  line  with  the  YZ  plane  to  its  specified  pitch  angle 
<p^,  meas\ired  between  the  pitch  chord  line  and  the  XY  plane.  Angles  measured 

counterclockwise  from  the  XY  plane  are  positive.  If  there  is  no  rake,  the  projection 
of  the  line  on  the  XZ  plane  runs  through  the  origin  of  the  XZ  axis  system  which 
becomes  the  center  of  pitch  rotation  (Figure  10). 

The  X  and  Z  coordinates  of  the  mid  chord  SV2  of  the  unwrapped  section 
(Figure  36a) 


ZK.  =  SX.  sin  4).  ♦  RK,  (2.35) 

111  i 

XK.  =  -  SK.  cos  4),  (2.36) 

1  1  i 

when  warp  is  specified  instead  of  skew  (Figure  14)  then 

ZK.  =  0 
1 

XKj  =  -  R^a^  =  -  WK^  (2.36a) 

The  shape  of  the  section  about  the  pitch  line  is  described  in  the  local 
coordinate  system  X’  Z*  (Figure  7)  by  the  thickness  offsets  TB^^  and  TF^j  at  fractions 

Cj  of  the  chord  line  (Figure  36a).  The  coordinates  of  the  face  and  back  can  be 

transformed  from  the  local  system  to  the  global  XYZ  system  using  the  following 
equations  and  remembering  that  TB^j  and  TF^j  are  negative  above  the  chord  line  and 

positive  below. 
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(2.37) 


where 


ZB..  =  ZK.  +  S,  (c.  -  Va)siii  4),  +  TB..  cos  d), 

1)  1  I  )  i  1) 


R.  =  RF^  -RR 


=  -XD../R. 
B  1]  1 

i) 


=  -  1  (XK  -  S  (c  Va)  cos  4>  .+  TB  sin  <t>  ) 

R7  i  i  j  1  i)  I 


(2.38) 

(2.39) 

(2.39a) 


Wrapping  the  section  about  the  specified  cylindrical  radius  R^  gives 


XB..  =  -  R^  sinog 


1) 


YB. .  =  R.  COSUn 

1)  1  B, 


(2.40) 

(2.41) 


RF. 


fraction  of  full  propeller  radius,  RR; 


c. 

) 


fraction  of  chord  length,  for  the  )th  offset  location; 


1 

TB. 


cylindrical  radius  for  section  i; 


back  of  blade  section  offset  at  Cj  of  radius,  R^; 


XB^.  =  X  coordinate  of  back  of  blade  section  offset,  TB^^; 

YB„  =  Y  coordinate  of  back  of  blade  section  offset,  TB^^; 

ZB..  =  Z  coordinate  of  back  of  blade  section  offset,  TB,.; 

1)  i) 


ZK. 


Z  coordinate  of  mid  chord; 


XK. 


X  coordinate  of  mid  chord  of  unwrapped  section; 


chord  length  of  section  at  R^; 


pitch  angle  of  chord,  S^,  at  R^,  positive  when  measured 
counterclockwise  from  the  X  axis; 


XDj,  a  projected  arc  length  to  TB^^  from  YZ  plane  (Figure  36b); 

(Also  X  coordinate  of  point  TB^^  of  unwraped  section  (Figure  36a)) 
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o  *  projected  angle  in  radians  to  offset  TB..  from  the  Y  axis  at  R.. 

®lj  ^ 

The  same  equations  are  used  for  the  face  (pressure  side)  coordinates  of 

the  blade  except  that  TF,,  are  substitutes  for  TB.,,  etc. 

ij  1) 

In  this  manner,  the  two  dimensional  data  supplied  by  the  designer  is 
converted  to  three  dimensional  cartesian  coordinates,  and  projections  of  the  blade 
on  the  XY,  XZ  and  YZ  planes  can  be  produced. 

2.13  Conventional  Views  of  the  Blade 

2.13.1  Transverse  View 

Once  the  propeller  blade  is  described  in  three  dimensional  cartesian 
coordinates,  various  views  can  be  produced  through  the  use  of  computer  graphics. 
One  conventional  view  is  a  projection  on  the  XY  plane  which  is  referred  to  here  as 
the  transverse  view.  It  is  a  view  of  the  blade  an  observer  would  see  if  he  were 
standing  astern  of  the  ship  looking  forward  at  the  propeller  set  at  its  design  pitch 
(Figure  37). 

2.13.2  Longitudinal  View 

This  is  a  projection  of  the  blade  on  the  YZ  plane.  It  is  a  view  an  observer 
would  see  when  looking  at  the  propeller  from  the  starboard  side  (Figure  38). 

2.13.3  Plan  View 

The  plan  view  is  not  always  given.  When  supplied,  it  is  a  projection  of  the 
blade  on  the  XZ  plane  (Figure  39). 

2.13.4  Expanded  View 

The  expanded  view  (Figure  11)  is  generally  Included  in  the  suite  of  graphic 
representations  of  the  propeller  geometry  as  it  provides  a  compact  method  of 
displaying  the  blade  sections  and  the  skew. 

2.13.5  Pitch  Diagram 

The  pitch  diagram  (Figure  9)  is  drawn  to  show  the  pitch  variation  with 
blade  radius.  It  should  form  a  smooth  curve  from  which  the  P/D  ratio  can  be 
determined. 

2.13.6  Blade  Thickness  Diagram 

The  blade  thickness  diagram  (Figure  40)  illustrates  the  blade  spanwise 
maximum  thickness  neglecting  rake  or  skew.  The  thickness  for  each  specified  blade 
section  is  shown  rotated  into  the  YZ  plane  to  give  the  true  scaled  dimension  for  the 
thickness  and  its  location  relative  to  the  Y  axes. 
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2.13.7  Combined  Blade  Thickness  and  Rake  Diagram 

The  applied  rake  of  the  blade  is  shown  in  Figure  41  in  combination  with 
the  maximum  thickness,  in  true  scaled  dimension,  on  the  YZ  plane.  Positive  rake  is 
aft  as  indicated  by  the  arrowhead  representing  the  propeller  shaft  axis.  The  angle 
formed  by  a  line  from  the  tip  to  the  point  of  intersection  of  the  shaft  axis  and  the  Y 
axis  is  the  global  rake  angle. 

2.13.8  Net  Rake  and  Blade  Thickness  Diagram 

The  net  rake  is  the  resultant  of  the  imposed  rake  plus  the  skew  induced 
rake.  It  is  shovm  in  Figure  42  in  combination  with  the  maximum  thickness,  in  true 
scaled  dimension,  on  the  YZ  plane. 

2.13.9  Combined  Drawring  of  Propeller  Blade 

The  conventional  views  of  the  blade  are  frequently  combined  to  provide 
all  the  blade  information  on  a  single  drawing  as  shovm  in  third  angle  projection  in 
Figure  43.  Thus,  the  transverse  view  superimposed  on  the  expanded  outline  and  the 
rake  and  thickness  superimposed  on  the  longitudinal  view  are  presented  with  the 
pitch  diagram.  Also,  the  outline  of  the  swept  area  of  interest  is  shown  superimposed 
on  the  longitudinal  view.  The  swept  area  is  especially  useful  information  to  have  for 
highly  skewed  blades  as  it  shows  the  area  that  the  blade  will  cut  from  the  YZ  plane. 
Required  clearances  around  the  blade  can  be  shown  in  this  manner. 

2.14  Rotation  of  Blade  Axes 


It  is  useful  to  have  the  ability  to  rotate  the  propeller  blade  and  geometric 
axes  system  to  new  positions  to  produce  a  variety  of  views  which  wUl  give  a  better 
representation  of  the  blade  than  straight  orthographic  projections.  Rotations  can  be 
about  either  the  Global  XYZ  axes  or  the  Local  body  xyz  axes  systems  using  the  same 
coordinate  axis  transformation  matrices.  However,  the  particular  axes  system 
selected  will  dictate  the  size  of  the  rotation  angles  and  the  order  in  which  the 
transformation  matrices  are  applied.  For  example,  in  accordance  with  the  right 
hand  rule,  isometric  drawings  can  be  produced  on  the  XY  plane  either  by  first 
rotating  -45°  about  the  vertical  global  Y  and  then  30°  about  global  X  or  by  rotating 
-45°  about  local  y,  -24.94°  about  local  z  and  24.94°  about  local  x.  In  the  latter  case, 
however,  the  local  axes  system  rotates  with  the  body  and  the  multiplications  of  the 
rotational  transformation  matrices  have  to  be  done  in  the  reverse  order;  first  the  x 
matrix,  then  the  z  matrix  and  finally  the  y  matrix.  A  discussion  of  the  three 
rotational  transformation  matrices  is  given  briefly. 


2.14.1  Rotation  4>  About  the  Global  Y  or  Local  y  Axis. 


The  global  coordinates  of  the  end  points  of  the  new  position  of  the 
original  xyz  propeller  body  axis  system  are  obtained  by  using  the  rotational 
transformation  matrix. 


[<t>]  * 


cos  d> 
0 

sin4> 


0 

1 

0 


-sin  4> 
0 


cos 


(2.42) 


Assuming  the  local  axes  are  unit  vectors,  the  new  global  coordinates  are  simply 
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({X)  (Y)  izj)  = 


l<t>] 


(2.43) 


The  same  transformation  matrix  can  be  applied  to  any  single  point.  With 
existing  xyz  coordinates,  the  new  coordinates  of  the  point  would  be 


X 

Y 

Z 


T 


[45] 


(2.44) 


2.14.2  Rotation  a  About  the  Global  Z  or  Local  z  axis 


is 


The  coordinate  transformation  matrix  for  rotation  about  Z  (or  z)  is 

(a)  *  - 

2.14.3  Rotation  6  About  the  Global  X  or  Local  x  Axis 

The  coordinate  transformation  matrix  for  rotation  about  the  X  (or  x)  axis 


cos  a  sin  a 
sin  a  cos  a 
0  0 


(2.45) 


10]  . 


1 

0 

0 


0  0 
cos  0  sinG 
sin  0  COS0 


(2.46) 


As  mentioned,  the  rotations  can  be  combined  and  the  resulting 
coordinates  of  the  geometry  referenced  to  the  global  axis  position.  The  sequence  of 
rotations,  however,  is  important  when  rotating  about  the  local  axes  because  the 
coordinates  are  obtained  by  working  back  from  the  final  axes  position.  Thus,  for 
rotations  about  the  y,  z  and  x  axes,  the  transformation  matrices  are  multiplied  in 
the  order  of  (0],  [a]  and  [4>].  The  net  mutiplication  [0]  [a]  [(t>)  of  the  three 
transformation  matrices  is  sometimes  referred  to  as  the  literal  matrix  [L]  =  (©]  (a] 
[4>]. 


3.  COORDINATES  FOR  MANUFACTURE  AND  MEASURKMKNT 

The  coordinates  that  are  usually  generated  are  for  the  propeller  blade  at 
the  design  pitch  angle.  The  pitch  angle  at  the  70%  radius  is  used  as  the  datum  for 
the  pitch  settings  when  a  variable  pitch  blade  is  rotated  about  its  spindle  axis.  At 
the  design  pitch  setting,  the  coordinates  can  be  used  for  stress  analysis  and  graphic 
representation  of  the  blade.  The  design  pitch  is  not,  however,  the  best  angle  for 
machining  with  a  numerically  controlled  milling  machine  or  for  measuring  the  blade 
surface  contours.  For  these  two  processes,  It  is  often  desirable  to  present  the 
maximum  projected  area  possible  to  reduce  machine  tool  and  measurement  gauge 
travel.  Zero  pitch  at  the  70%  radius  gives  close  to  the  maximum  projected  area  for 
most  blades.  Other  angles  may  be  used  if  they  are  found  to  be  more  suitable. 
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3.1  Cartesian  Coordinates 


To  present  the  face  of  the  blade  for  measurement  or  machining,  it  is 
rotated  about  its  spindle  axis  through  the  angle  -  <>  to  approximately  zero  pitch. 
Right  handed  blades,  i.e.,  blades  which  advance  when  the  propeller  shaft  is  rotated 
clockwise,  are  rotated  clockwise  about  their  spindle  axis.  Left  hand  blades  are 
rotated  counterclockwise.  The  rotation  -  4>  reqiaires  a  coordinate  transformation  of 
the  same  form  as  was  described  in  Section  2.14.1.  The  values  of  YB  and  YF  are  not 
changed  by  the  rotation.  The  values  of  XF,  XB  and  ZF  and  ZB,  however,  are 
modified  as  follows 


XMF,.  =  XF,,  cosd)  4  ZF,,  sin4> 
i)  ij  ^  i) 

(3.1) 

ZMF^.  =  “XF^j  sin0  4  ZF^^  cos4> 

(3.2) 

To  gain  acces:  to  the  back  surface,  the  blade  must  be  rotated  an 
additional  180®.  Thus, 


XMB  .  =  XB,  cos(4)+tt)  ZB,,  sin(d>+TT)  (3.3) 

i)  i)  i) 

ZMB..  =  -XB.,  sin(<t>+Tr)  4  ZB,,  cos(<b4Tr)  (3.4) 

i)  i)  i) 

The  coordinates  of  the  blade  are  now  in  a  form  that  permits  easy 
measurement  of  the  blade  front  and  back  surfaces  provided  the  measurement  device 
can  measure  accurately  along  three  axes,  X,  Y  and  Z.  An  N.C.  milling  machine  or  a 
}ig  borer  are  devices  that  can  be  \ised  for  this  purpose.  A  typical  example  of  blade 
geometry  in  cartesian  coordinates  form  is  shown  in  Table  1. 

3.2  Cylindrical  Coordinates 


The  blade  to  this  point  has  been  described  in  spatial  cartesian 
coordinates.  Cylindrical  coordinates  are  better  suited  for  measurement  by  a  gauge 
of  the  type  shown  diagramatically  in  Figure  44.  For  this  purpose,  the  coordinates 
are  converted  to  an  angle  measured  from  the  Y  axis,  a  radius  R  from  the  center  of 
rotation  of  the  shaft  of  the  propeller  and  a  Z  axis  dimension  for  every  section  offset 
given  in  the  designer's  original  data.  The  measurement  angles  are 

Face  Angle 

=  -  tan'^  (XMF,,  /YF,,)  (3.5) 

F  1)  i} 

Back  Angle 

Bg  .  -  tan'^  (XMB^j/YBjj)  (3.6) 
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Face  Radii 


RMF,,  =  >/XMF,.^  +  YF,,^  (3.7) 

ij  V 

Back  Radii 


RMB^j  =  +  YBjj*  (3.8) 

The  measurement  of  the  Z  coordinates  can  be  more  complex  if  a  dial 
gauge  is  used  and  negative  values  cannot  be  read  directly.  If  negative  values  occur 
in  the  column  of  Z  offsets,  an  origin  shift  is  required  to  convert  all  the  offsets  to 
positive  values.  This  is  accomplished  by  subtracting  the  largest  negative  Z  from  all 
the  others  (Figures  44  and  45).  The  subroutine  RXYZ  of  program  BLADGM  carries 
out  the  origin  shift  when  required  and  prints  out  the  modified  Z  offsets  as  a  data 
column  titled  "GAGE". 


4.  COMPUTER  PROGRAMS 


A  computer  program,  BLADE,  consisting  of  three  subprograms  BLADAT, 
BLADGM  and  BLADPL  has  been  developed  for  generating  and  displaying  propeller 
blade  geometry. 

BLADAT  is  an  interactive  program  which  prompts  the  user  in  entering 
propeller  geometry  data.  It  can  accept  the  data  in  the  three  following  basic  forms 
which  are  used  to  describe  the  blade  sections. 

a.  blade  section  data  in  non-dimensional  form  as  ratios  of  section  half 
thickness  over  maximum  thickness  and  camber  over  maximum  camber  at 
specified  fractions  of  the  section  chord. 

b.  NACA  section  data  in  the  form  of  basic  and  mean  line  data. 

c.  section  offsets  from  drawings  or  tables. 

Upon  completion  of  the  basic  data  entry,  the  data  can  be  examined  by 
plotting  or  by  directly  examining  the  input.  At  this  stage  any  errors  which  are 
identified  can  be  corrected  using  BLADAT. 

BLADGM  is  also  an  interactive  program.  It  accepts  the  data  from 
BLADAT  and  completes  the  blade  section  description  by  generating  the  additional 
leading  and  trailing  edge  coordinates  when  required.  It  also  completes  the  data 
required  to  describe  the  blade  tip  and  generates  the  three-dimensional  coordinates 
,in  tabular  form  for  the  propeller,  required  for  blade  manufacture  and  checking. 
The  three-dimensional  coordinates  in  this  form  can  also  be  used  to  generate  finite 
element  models  of  the  blade  for  strength  and  vibration  analysis. 
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The  subprogram  BLADPL  \jses  the  data  generated  by  6LADGM  to  plot  the 
propeller  geometry.  It  provides  a  menu  of  the  plotting  options  for  selecting  the 
graphics  required  to  produce  conventional  dravdng  of  an  individual  blade  or  a 
complete  propeller.  In  addition,  it  offers  an  option  to  view  the  blade  in  a  solid 
geometry  form  by  generating  a  session  file  for  the  solid  modelling  program 
PATRAN.  An  example  of  views  of  blade  geometry  generated  by  PATRAN  from  a 
BLADPL  produced  session  file  is  shown  in  Figure  46. 

The  versatility  of  the  program  is  shown  by  its  ability  to  produce  drawings 
of  a  variety  of  blade  forms.  Figure  47  is  the  plan  view  of  the  blade  sections  of  a 
moderately  skewed  blade.  Figure  48  is  the  same  view  of  a  long  narrow  blade  with 
high  skew  and  Figxire  49  shows  a  transverse  projection  of  a  six  bladed  propeller  with 
moderate  skew. 


5.  CONCLUSIONS 


The  geometry  of  marine  propellers  has  been  described  and  the 
relationships  have  been  derived  for  generating  propeller  blade  geometry  in  the  form 
of  three  dimensional  coordinates  from  two  dimensional  blade  section  data  including 
pitch,  rake  and  skew. 

The  problem  of  curve  fitting  the  leading  and  trailing  edges  and  the  blade 
tip  has  been  overcome  by  use  of  a  special  coordinate  transformation  and  cubic  spline 
interpolation  for  non  elliptic  and  elliptic  form. 

The  equations  derived  for  defining  the  geometry  of  a  blade  have  been  used 
to  develop  a  computer  program  for  blade  coordinate  generation.  The  program 
inclxodes  the  ability  to  produce  computer  generated  graphic  displays  of  blade 
geometry  for  checking  hydrodynamic  form  and  the  detection  of  data  errors.  It  also 
can  be  used  for  producing  conventional  drawings  and  tabular  data  for  manufacture 
and  blade  form  measurement. 
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TABLE  1 :  COORDINATES  OF  THE  BLADE  BACK  SURFACE 
ROTATED  FOR  EASY  MEASUREMENT 


3  DIMENSIONAL  X.Y.Z  CARTESIAN  COORDINATES  FOR  CHECKING. 

BLADE  ROTATED  180  DEGREES  FROM  FACE  SETTING  OF  0.000  DEGREES 
FRACTION  RADIUS  0.3000  BLADE  RADIUS  502.95  MM. 


FRACTION 

BACK 

CHORD 

ZMB 

YB 

ZMB 

0.0000 

-188.13 

466.44 

363.87 

0.0010 

-185.16 

467.63 

364.52 

0.0038 

-181.34 

469.12 

363.94 

0.0075 

-177.01 

470.77 

362.67 

0.0100 

-174.35 

471.76 

361.70 

0.0250 

-161.34 

476.37 

354.24 

0.0500 

-143.07 

482.17 

339.89 

0.1000 

-110.44 

490.67 

308.87 

0.1500 

-81.31 

496.33 

275.82 

0.2000 

-55.09 

499.92 

241 . 16 

0.2500 

-30.27 

502.04 

205.72 

0.3000 

-6.93 

502.90 

169.49 

0.3500 

15.01 

502.73 

132 . 56 

0.4000 

35.70 

501.68 

95.01 

0.4500 

55.20 

499.91 

56.89 

0.5000 

73.54 

497.54 

18.23 

0.5500 

90.66 

494.71 

-21.00 

0.6000 

106.46 

491 . 55 

-60.87 

0.6500 

120.93 

488.20 

-101.35 

0.7000 

134.15 

484.73 

-142.44 

0.7500 

146.12 

461.26 

-184.13 

0.8000 

156.84 

477.87 

-226.43 

0.8500 

166.44 

474.61 

-269.29 

0.9000 

174.92 

471.65 

-312.71 

0.9500 

182.60 

468.63 

-356.53 

0.9750 

186.08 

467.26 

-378.62 

0.9875 

187.39 

466.74 

-389.90 

0.9967 

187.96 

466.51 

-398.37 

1.0000 

187.01 

466.89 

-402.10 
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FAN  CIRCULAR  ELLIPTIC 


WEEDPROOF  WIDE  SKEW-BACK  OVOID  OVOID  DAGGER 


FIGURE  1 :  PROFILE  SHAPES  FOR  SCREW  PROPELLERS  (REF.  1 ) 


HGURE  2:  PARTS  OF  A  MARINE  PROPELLER 
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CYLINDRICAL 


SINGLE  RADIUS  FILLET  COM 


»lim» 
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FIGURE  6:  PROPELLER  AXES  SYSTEMS 


FIGURE  7:  FULLY  DEVELOPED  SECTION  PROJECTED  ON  XZ 
PLANE  WITH  ORIGIN  SHIFTED  TO  MIDCHORD 
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FIGURE  8; 


PITCH  OF  AN  UNWRAPPED  BLADE  SECTION  BY  DEVELOPMENT 
OF  AN  XZ  PLANE  FROM  CYLINDRICAL  SURFACE 


nCURE  9:  PITCH  DIAGRAM 
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SKi  COS^ 


FIGURE  10:  PROPELLER  SECTION  UNWRAPPED  ON  THE  XZ  PLANE 


HGUREll:  EXPANDED  VIEW  OF  BLADE 


-31  - 


FIGURE  12:  WRAPPED  BLADE  SECTION  SHOWING  SKEW  ANGLE  ai 


FIGURE  13:  FULLY  DEVELOPED  SECTION  ON  XZ  PLANE 
SHOWING  SKEW  AND  RAKE 


a)  UNWRAPPED  VIEW  OF  SECTION  SHOWING 

RELATIONSHIP  OF  WARP  AND  SKEW 


RGURE  14:  BLADE  WARP 
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b)  WARP  ANGLE  SHOWN  ON  PROJECTED  VIEW  OF  BLADE 


c)  LINEAR  DIMENSION  FOR  WARP  AND  SKEW 
ON  EXPANDED  MIDCHORD  DATUM 


HGUREU:  BLADE  WARP 
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a)  EXPANDED  BLADE  SECTION 


FIGURE  1 5 :  INITIAL  DESCRIPTION  OF  BLADE  SECTION 


FACE  PITCH  LINE 


FIGURE  1 6;  FACE  PITCH  LINE  FOR  FLAT  OR  CONCAVE  SECTIONS 
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CHORD  LINE 


FIGURE  17:  CONVERSION  OF  MEAN  LINE  AND  BASIC  FORM 
DATA  TO  SECTION  OFFSETS 


FIGURE  18:  LEADING  EDGE  GEOMETRY 


-36- 


CURVE  FITTING  INTERVALS  FOR  LEADING  EDGE  HT 


b)  TRANSFORMATION  OF  X-AXIS  TO  ^-AXIS 


FIGURE  20:  ELLIPTICAL  LEADING  EDGE 


FIGURE  2 1 :  LEADING  EDGE  RADIUS  FITTED  TO  CAMBERED  SECTION  DATA 
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b)  CHORD  STATIONS  SPECIFIED  AS  FRACTIONS 

FROM  POINT  OF  MAXIMUM  THICKNESS 


HGURE  22;  TRANSFORMATION  OF  SECTION  COORDINATES 
USING  CHORD  FRACTIONS 
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FIGURE  23:  CURVE  FITTING  BLADE  TIP 
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FIGURE  24:  BLADE  SPANWISE  MAXIMUM  THICKNESS  SECTION 


-40- 


EXPANDED  SECTIONS 


FIGURE  25:  ADDITIONAL  BLADE  SECTIONS  GENERATED  FROM 
LAST  GIVEN  SECTION  TO  THE  TIP  OF  THE  BLADE 


nGURE26:  COMPOUND  RADIUS  TYPE  FILLET 
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BUDE  SECTION  RADIUS 


Inches 


nOURE  27:  DATA  PLOT  OF  PITCH  ANGLE  VERSUS  BLADE  SECTION  RADIUS 
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Inches 


FIGURE  28:  DATA  PLOT  OF  SKEW  VERSUS  BLADE  SECTION  RADIUS 
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BLADE  SECTION  RADIUS 
Mi  1 1 i««t«rs 


FIGURE  3 1 :  DATA  PLOT  OF  LEADING  EDGE  RADIUS 
VERSUS  BLADE  SECTION  RADIUS 
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HGURE  32:  DATA  PLOTS  OF  TRAILING  EDGE  RADIUS 
VERSUS  BLADE  SECTION  RADIUS 
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BLADE  SECTION  RADIUS 
Inches 


FIGURE  34:  DATA  PLOT  OF  MAXIMUM  THICKNESS 
VERSUS  BLADE  SECTION  RADIUS 
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FIGURE  35:  INDIVIDUAL  SECTION  PLOT  WITH  BACK  AND  FACE  OFFSETS 
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TRANSVERSE  PROJECTION 


FIGURE  37:  COMPUTER  PLOT  OF  TRANSVERSE  VIEW 


LONGITUDINAL  PROJECTION 


HGURE  38:  COMPUTER  PLOT  OF  LONGITUDINAL  VIEW 


MAXIMUM  THICKTCSS  WITHOUT  RAKE  OR  SKEW 


FIGURE  40:  BLADE  THICKNESS  WITHOUT  RAKE  OR  SKEW 
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HGURE  41 :  COMBINED  BLADE  TfflCKNESS  AND  RAKE 
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FIGURE  47:  PLAN  VIEW  OF  ASSEMBLED  BLADE  SECTIONS 
FOR  A  MODERATELY  SKEWED  BLADE 
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HGURE  48:  PLAN  VIEW  OF  ASSEMBLED  BLADE  SECTIONS 
FOR  A  HIGHLY  SKEWED  BLADE 


-62- 


HGURE  49: 


TEIANS  VERSE  VIEW  OF  THE  COMPLETE 
MODERATELY  SKEWED  PROPELLER 
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APPENDIX  A 


TRANSFORMATION  FOR  LEADING  EDGE 


TrAnsformation  of  the  x  axis  in  the  interval  x  ^  x  ^  x  using  x  s 

1  2 

wrill  overcome  the  problem  of  infinite  slope  at  x^  and.  therefore,  infinite  value  of 
the  slope  coefficient  for  the  x  cubic  spline  function. 

The  transformation  x  »  'was  arrived  at  because  of  the  approximation 
of  the  nose  of  the  blade  section  to  the  arc  of  a  circle  with  leading  edge  radius  r  (see 
also  Figure  19a  and  19bl. 

The  equation  for  a  circle  with  its  centre  on  the  x  axis  at  radius  r  is 

(x-r)*  +  y*  =  r’  (A.l) 

Expanded  and  rearranging 

y^  =  2rx  -  x^  (A.la) 

for  small  x  values  (ie.  equivalent  to  nose  of  blade  section) 

y^~2rx  (A. lb) 

or  taking  )ust  the  positive  root 

y  ~  VTr  vTC  (A.lc) 

Letting  -/x  =  C  (ie-  transformation  x  =  C*)  *  linear  relationship  is  produced. 

y  *  VTt  f  (A.2) 

Resulting  in  finite  slope  at  C  -  0  of 

42  =  vTr  (A.2a) 

dC 

In  the  interval  ^  ^  (ie.  equivalent  to  x^  $  x  $  x^)  a  cubic  spline 

function  is  used  for  interpolation  and  to  generate  additional  ordinate  values. 

y  =  A  ♦  B  (C-C  )  +  C  (C-C  )’  ♦  D  (C-C  )*  (A.3) 

1  X  1  X  1  X  X 


-  65- 


where,  at  C  =  =  0 

y  =  A  =0 
1 

dX  =  B  *  VTF 
dC 

and  at  C  =  C 

2 

y  =  B  C  *  C  C’  +  D  C*  (A.4a) 

12  12  12 

4x  =  B  +  2C  f  +  3D  (A.4b) 

df  1  12  12 

by  equating  the  offset  ordinate  and  slope  given  by  the  last  two  equations  to  the 
values  for  the  next  interval  (see  Figure  19(a)  and  19(b))  such  as 

y  =  A^ 

^  ^  ^  =  2B  C  (A.5) 

dC  dx  dC  *  * 

and  using  the  transformation  C  =  » the  coefficients  A  ,  B  ,  C  and  D  can  be 

X  1  X  1 

obtained 


=  0 

(A.6) 

=  V7r 

(A.7) 

,  1 

II 

==A  -  aB  vir  -  2B 

(A.8) 

=  JL 

X 

2 

[Ib  VTT  -  ♦  bI 

(A.9) 
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APPENDIX  B 


CURVE  PITTING  THE  BLADE  TIP 


To  determine  the  radius  of  curvature  at  the  tip  of  the  blade,  a  parabola  is 
fitted  of  the  form  (Figure  23) 

y  =  ax*  (B.l) 

The  radius  of  curvature  is  equal  to  the  reciprocal  of  the  curvature  K  of  a  general 
curve  y  =  f(x)  given  by 


r  =-L=  [  ..UCdy/dx)*]  (B.2) 

*=  K 

Evaluating  the  first  and  second  derivatives  of  y  at  x  =  0,  the  radius  of  curvature  is 
given  by 

r^  =  i  (B.2a) 

c  2a 


Using  the  coordinates  x  and  y  of  the  leading  and  trailing  edge  points  of  the  last 
specified  blade  section, 


a  s 


(B.2b) 


Thus,  the  radius  of  the  blade's  tip  is 


(B.3) 


This  radius  is  used  to  determine  the  slope  coefficient  of  the  cubic  spline  function  in 
the  transformed  coordinate  system  for  the  blade  tip  region  similar  to  that  outlined 
in  Section  2. 7. 3.1  and  Appendix  A  for  the  leading  edge  fit. 


-67- 


REFERENCES 

1.  Sa-unders,  Harold,  Hydrodynamics  in  Ship  Design,  Vol.  1,  SNAME  461-2,  1957. 

2.  O'Brien,  T.P.,  The  Design  of  Marine  Screw  Propellers,  Hutchinson  &  Co.,  Ltd., 
London,  1962. 

3.  Abbott,  H.I.  and  Von  Doenhoff,  A.E.  Theory  of  Wing  Sections.  Dover 
Publications,  Inc.,  New  York,  1959. 

4.  Gerald  Curtis,  F.,  Applied  Numerical  Analysis,  Addison- Wesley  Publishing  Co. 
Massachusetts. 


-68- 


UNCLASSIFIED 


SECURITY  CLASSFCATION  OF  FORM 
(highMt  daitification  of  Titi*.  Abatraet.  K*y«»ord$) 


DOCUMENT  CONTROL  DATA 

(Security  classHicatior^  of  titlo.  body  of  abstract  artd  iodaiing  anocMation  must  ba  antorad  wban  tha  ovarall  documant  is  ciassifiad) 


1  ORIGINATOR  (tha  narrw  and  address  of  Uta  organtzaiion  prapshng  tha  document 
Organizations  for  whom  the  document  was  prepared,  a.g.  Establishmant  spenaoring  a 
contractor’s  report,  or  tasking  agency,  are  antarad  in  aaetion  A.) 

2.  SECURITY  CLASSIFICATION  (ovwall  Mcufty 
dauifioation  el  tha  documant,  including 
ipacial  warning  tarmc  H  appiicabla) 

Defence  Research  Establishment  Atlantic 

UNCLASSIFIED 

3  TITLE  (tha  corrpiata  docurrwnt  titia  as  indieatad  on  tha  titia  p^a.  lU  otassilioation  should  ba  indicatad  by  tha  apprepriala  abbra^ion 
(S.C.R  or  U)  in  paranthasas  after  tha  title.) 

THE  GEOMETRY  OF  MARINE  PROPELLERS 


4.  AUTHORS  (Last  name,  first  name,  middle  initial.  K  milttary.  show  rank.  a.9.  Doe.  Maj.  John  E.) 

Smith,  Donald  R.  and  Slater,  John  E. 


5  DATE  OF  Publication  (month  and  year  of  publicaton  of  doeumanl}  PAGES  (total  oontaimrtg  NO.  OF  REFS  (total  oitad  in 

information.  Inoluda  Annaxas,  documar^t) 

November  1988  Ap»»Bdi<».,  «c.)  4 


g  DESCRIPTIVE  NOTES  (the  category  of  the  documant.  a.g.  technical  report,  technical  note  or  memorandum.  H  appropriala.  enter  tha  type  of 
report,  a.g.  irttarim,  progress,  summary,  annual  or  final.  Oiva  tha  indushra  da.«>s  arhan  a  spacHie  tiporlirtg  period  is  eowarad.) 

Technical  Memorandiim 


3  SPONSORING  ACTIVITY  (the  name  of  tha  department  pro^  oHioa  or  laboratory  sponsoring  tha  raaaarch  and  davalopmartt.  Induda  tha 
address.) 

Defence  Research  Establishment  Atlantic 
Dartmouth,  Nova  Scotia 


PROJECT  OR  GRANT  NO.  (d  appropriate,  tha  applicabta  research  and  ^  CONTRACT  NO.  ( if  appropriate,  tha  appl^abia  number  under 
davaiopmant  project  or  grant  number  under  which  the  documant  was  which  tha  documant  was  written) 
written.  Please  specify  whether  project  or  grant) 


10a  ORIGINATOR  'S  DOCUMENT  NUMBER  (tha  officia/  documnt 

number  by  which  tha  documant  is  identified  by  tha  originoting 
activity.  This  number  must  ba  uniqud  to  this  document.) 

DREA  Technical  Memorandum  88/214 


10b.  OTHER  DOCUMENT  NOS.  (any  other  nuirbars  which  may  ba 
assigned  this  document  akhar  by  tha  ohginator  or  by  tha 
sponsor) 


11.  DOCUMENT  AVAILABILITY  (any  limitalions  on  furthor  dissomination  of  tha  documant  olhar  than  thoea  imposad  by  sacurity 
clsssification) 

)  UnlimHad  distribution 

(  )  Distribution  Kfrstad  to  dafanoa  departments  and  dafanoa  oontraotors;  further  dbtnbution  only  as  approved 

(  )  Distribution  limHod  to  dafanoa  dapartrr>anti  and  Canadiai  dafanoa  oontradors;  furthar  dittrfeutien  only  asipprevod 

(  )  Distribution  limitad  to  govornmant  daparlmants  and  agancias;  furlhar  distribution  only  as  apprevad 

(  )  Distribution  limHad  to  dafanoa  departments;  furthar  distrbMion  only  as  approved 

(  )  Other  (plaasa  spady): 


12..  DOCUMENT  ANNOUNCEMENT  (any  limitations  to  tha  bk>fiegraphic  symounoamant  of  this  doeumanl  This  wifi  normally  oorrasporwf  to  tha 
Oocumar^t  Availability  (11).  Howavar.  whara  futhar  distribution  (beyond  the  audianoa  apadfiad  In  11)  is  posskria,  a  widar  announoamant 

.ud-nc  m.y  b.  UNLIMITED 


-69- 


DC003  2m/87 


KCURITV  CLAStiriCATIOM  Of  ro«M 

13.  ABSTRACT  I  •  brief  aid  fictiitl  sumnnry  of  Ifit  4ee<«nefll  It  may  elM  ipbea  ilsewbere  ifl  ibe  My  el  the  becumeni  itself  It  is  highly 
desireble  that  the  ^street  of  classified  dacuments  be  unelassifiad  Each  garagrM  af  the  abstract  shall  begin  with  an  indication  of  the 
security  classification  of  the  mformation  ui  the  paragraph  lunless  the  document  itself  is  unclassified)  represented  as  (S).  (Cl.  Oil.  or  (Ul 
It  IS  nat  necessary  to  rnchide  here  abstracts  m  both  affical  languages  unless  the  test  is  bilinguall 


The  geometry  of  marine  propellers  is  described  in  this  report.  The 
types  of  data  used  for  blade  section  description  are  covered  and  methods 
for  generating  blade  section  leading  and  trailing  edge  forms  are  derived. 
A  method  for  generating  the  blade  tip  from  minimal  data  is  given.  The 
computer  program  BLADE  developed  for  generating  3  dimensional  propeller 
blade  coordinates  is  described  together  with  examples  of  its  graphic 
output . 


14  KEYWORDS  DESCRIPTORS  O'  IDENTIFIERS  ItechnicaHy  meaningful  terms  er  shen  phrases  that  charactenie  a  dccumem  am)  could  tt 
heloful  in  cataloguing  the  document  They  should  be  selected  so  that  no  security  classification  is  required  Identifiers,  such  as  equipmem 
model  designation,  iredt  name,  mihury  protect  code  name,  geographic  location  may  also  be  included  If  possible  keywords  should  be  selected 
from  a  published  thesaurus  eg.  Thesaurus  of  Engineering  and  Scientific  Terms  (TEST)  wid  that  thesaurus-identified  H  it  is  not  possible  to 
select  indening  terms  which  re  Unclassified,  the  classification  of  each  should  be  indicated  as  with  the  title.) 

Protiellers 

Marine 

Geometry 


UNCLASSIFIED 

MCURITV  CLASSIklCATION  Of  fObka 


-70- 


